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ABSTRACT 



Cosmic microwave background (CMB) polarimetry has the potential to provide revolutionary advances in cosmology. 
Future experiments to detect the very weak B mode signal in CMB polarization maps will require unprecedented 
sensitivity and control of systematic errors. Bolometric interferometry may provide a way to achieve these goals. In a 
bolometric interferometer (or other adding interferometer), phase shift sequences are applied to the inputs in order to 
recover the visibilities. Noise is minimized when the phase shift sequences corresponding to all visibilities are orthogonal. 
We present a systematic method for finding sequences that produce this orthogonality, approximately minimizing both 
the length of the time sequence and the number of discrete phase shift values required. When some baselines are 
geometrically equivalent, we can choose sequences that read out those baselines simultaneously, which has been shown 
to improve signal to noise ratio. 
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i_i 1. Introduction 



■ The field of observational cosmology has been advancing quickly in recent years. Observations of the cosmic microwave 
J> ■ background (CMB) radiation ha ve been leading the way, as evidenced by WMAP's highly successful mapping of CMB 

C<~) | anisotropy (Hin shaw et al.|[2008l ). DASI's detection of the polarized component of the CMB (|Leitch et al.ll2005f ), and the 

■ 2006 Nobel Prize in Physics awarded to Joh n Mather and Ge orge Smoot. Momentum is building for experiments that 
' characterize the CMB polarization in detail (jBock et al.ll2006h . 

(N ■ A linear polarization map can always be expres sed as the sum of two component maps, denoted E and B 
OO ' (|Seliak fc Zaldarriagalll997t iKamionkowski et al.lll997l ). CMB experiments to date have detected only the "curl- free" 
E component, which is produced primarily by (scalar) density perturbations. The "divergence- free" B component is not 
00 , produced by scalar perturbations at linear order, and is therefore a clean probe of other, smaller effects. In particular, 

■ inflationary models predict a B-mode signal produced by gravitational wave (tensor) perturbations in the early universe. 
| These B modes promise to hold key information about the process of inflation and particle physics above the Grand 

*h . Unification scale. The challenge of finding the B modes is no small task, however: the B component is expected to be at 
)\ ' least an order of magnitude weaker than the E component (which is itself small compared to the temperature anisotropy) 
, over all angular scales. Experiments to search for B modes will require unprecedented sensitivity and control of systematic 

- - i errors. 

Bolometric interferometry is one proposed method for a chieving these goals ()Polenta et al.l l2007t ICharlassierl 120081 : 
iTimbie et al.l 2006; K orotkov et al.ll2006t iTucker et ai1l2008T ). A bolometric interferometer is a marriage between highly 
sensitive, inc oherent bolomet r ic det ectors and the phase-sensitive, systematic-error-reducing observing technique of in- 
terferometry. Hamilto n et al.l (|2008f) have shown that a bolometric interferometer can achieve sensitivities comparable 
to traditional technologies. The question of whether bolometric interferometry is useful for CMB polarimetry will thus 
depend on the method's ability to cont rol systema tic errors. Systematic errors in interferometers are certainly differ- 
ent from those in imaging experiments (jBunnl 120071 ): it can be argued that interferometers are superior in this regard, 
although this question requires further research. 

In a bolometric interferometer, the signals from a set of input fecdhorns are combined with either a Butler combiner 
or a quasi-optical (Fizeau) combiner. In either case, bolometers measure the total power in the combined beam - that 
is, each bolometer is illuminated by signals from all of the input horns. Since the signal in each detector is the sum of 
all the inputs, a bolometric interferometer is an example of an "adding" interferometer (as opposed to traditional radio 
interferometers, which are "multiplying" interferometers) . One of the keys to making this method work is to arrange for 
the phase information to be encoded in the bolometer signals, so that individual pairwise visibilities can be extracted. To 
achieve this goal, a sequence of phase shifts can be applied to each of the input horns, in such a way that each visibility is 
phase-shifted in an independent fashion. The resulting time series can be solved for the individual visibilities. The phase 
shift sequences should be chosen so that this inversion can be done with minimal noise. 
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We would like the length of the phase shift sequence to be as short as possible, to avoid error due to 1// noise in the 
detectors. Clearly the number of phase shifts must be at least as large as the number of visibilities to be recovered. In 
the most general case, an TV-horn interferometer has N(N — l)/2 distinct visibilities, requiring long phase shift sequences 
for interferometers with many inputs. On the other hand, if the input horns are arranged in a regular pattern, such as 
a square array, then many antenna pairs correspond to identical visibilities. These can be given identical phase shifts 
and read out together. This coherent treatment of equivalent (or redundant) baselines has two advantages. First, it 
allows for shorter phas e shift sequences. Second, by coadding equivalent signals, the signal-to-noise ratio is improved 
(jCharlassier et al.ll2008l . hereinafter C08). 

C08 gave an excellent overview of how a bolometric interferometer works and considered the choice of phase shift 
sequences in detail. For the case of a square array of horns, the paper presented a method of phase modulating the inputs 
that gives equivalent baselines identical phase shift sequences. In this paper we present independently-developed work on 
phase modulation and coherent addition of baselines that complements the methods of C08. We consider general horn 
arrangements as well as a regular square lattice. In the general case, we present a method for finding the optimal phase 
shift sequence assuming no baselines are redundant. In the case of a square array, we present a refinement of the method 
of C08. Unlike the original method, which achieves optimal noise performance only in the limit as the number of time 
steps tends to infinity, our method is optimal for sequences of nearly or exactly the minimum possible length. 

In Sect. [2] we present our formalism for denoting sequences of phase shifts and consider the criterion for an optimal 
phase shift sequence. Section [3] introduces a shorthand notation and applies this to a method for constructing bases for 
phase sequences and selecting optimal sequences for interferometers without equivalent baselines. In Sect. [H we consider 
a square array of horns, accounting for redundant baselines. For simplicity, we consider only one linear polarization state 
in these sections, but in Sect. [5] we generalize to an array with two polarizations. Section [6] contains a brief concluding 
discussion, and Appendix A contains a useful mathematical result. 



2. Formalism 

Suppose that our interferometer has N input horns, each of which receives one electric field component. (We will consider 
the case where both x and y components are received in Sect. [5]) Assuming monochromatic radiation of frequency v 
for simplicity, the signal entering the jth horn can be written Eje vt . This signal is presumed to have already been 
averaged over the antenna pattern, so no integral over the sky position will be explicitly shown in the equations below. 
We apply a time-dependent phase shift <pj (t) to each of these inputs. Since bolometers measure the total power, the 
signal detected by the mth detector is 



S m (t) oc 



JV 



J2 E E* k e l ^-^ + ^ {t) ~* k{t) \ (1) 

3,k=l 



where the phase shifts A jm are fixed by the geometry of the system and do not vary in time. In this expression, the 
detector can correspond either to a single point on the focal plane of a quasioptical combiner or to a single output of a 
Butler combiner. 

We are assuming here that all inputs contribute to all detectors with equal amplitude. If this assumption is relaxed, 
then an additional real factor Aj m would need to be included in each term of the sum. The presence of these factors 
would affect the overall sensitivity of the detector to the various visibilities, but we do not expect it to influence the 
optimal choice of phase shifts, so we omit it . 

Let us assume for the moment that we wish to recover the visibility associated with each pair of horns separately; we 
will return below to the case in which redundant baselines are coherently added before detection. We wish to choose the 
phase shifts 4>j(t) to enable recovery of all of the cross terms EjE% from each detector. In principle, we could aim for a 
weaker goal, namely to ensure that each cross term be recoverable from the full set of detector outputs Si(t), S^i), . . .; 
however, to avoid systematic errors resulting from subtracting signals in different detectors, it is preferable to insist that 
each visibility be recovered from each detector separately. 

Since we are focusing on one detector at a time, we suppress the subscript m in eq. ([T]). Furthermore, the time- 
independent phase shifts A jm do not affect the problem of visibility recovery, so we suppress these as well. Finally, we 
assume that the phase shifts <pj are changed in discrete time steps, so we replace the functions <f>j(t) with sequences <f>j t . 
Here t = 1, 2, . . . , M, where M is the number of steps in the phase shift sequence at each detector. With these changes, 
equation ([1]) becomes 

N 

S t oc EjEijfa*-*") = Vjke*^*-^. (2) 

j,k=l 

Finally, we assume that the phase shifts can take on P equally-spaced values from to 2ir: 

<f> jt € {2irp/P | p = 0,l,2,...,P-l}. (3) 

Given the time sequence of measurements Si, S2, ■ ■ ■ , Sm, we wish to recover all N(N — 1) complex visibilities 
Vjk = EjE*k w ith j 7^ fc (or equivalently to recover both real and imaginary parts of all pairs with j < k). In addition, 
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we will always recover the total power term J^. Vjj = • l-E^I 2 , which enters each St equally. Solving for the cross terms 
is therefore simply inverting a linear system of M equations for N(N — 1) + 1 unknowns. Generically, we expect this to 
be possible as long as 

M > N(N - 1) + 1 = N vis . (4) 

We want to insist not just that the visibilities be recovered, but that they be recovered with minimal possible noise. 
To be specific, the recovery problem we wish to solve is 

S = AV, (5) 

where S is the M-dimensional signal vector, V is the iV v i s -dimensional vector of visibilities to be recovered, and A is a 
matrix whose elements are determined by the phase shifts: 

A tm = e i ^*-***\ (6) 

where the mth visibility V m corresponds to the horn pair jk. In this situation, where all elements of the matrix A have 
absolute value equal to one, the minimum possible noise contributions to the visibilities is achieved when all columns of 
A are orthogonal: 



A* tm A tm ' = for m + ml (7) 



or equivalently for A^A proportional to the identity matrix. 

A proof of this statement is provided in Appendix |XJ Intuitively, it says that the visibilities are recovered with 
minimum noise when they are maximally independent of each other, that is, when they contribute orthogonally to the 
time series of signals at the detector. 

Let us summarize. For any pair of horns j, k, define an M-dimensional vector 

Our goal is to choose the set of phase shifts 4>jt such that the vectors <frjfc and Qyk 1 are orthogonal whenever (jk) ^ (j'k 1 ). 
When this condition is satisfied, each visibility is recovered simply by taking the dot product of the detector signal with 
the corresponding vector <&: the estimator of Vjk is 

*» = M*» 8= M 

t=i 

We will call the vector the "mask" for the baseline jk. 

Note that 4"^ and &kj are complex conjugates of each other. The requirement that these be orthogonal, which means 
roughly that the elements of uniformly sample directions in the complex plane, is necessary for both the real and 
imaginary parts of Vjk to be recovered with minimum noise. 

It may be instructive to compare the phase shift schemes for the bolometric interferometer with those applied in 
a traditional multiplying interferometer. In traditional interferometry, orthogonal patterns of square-wave phase shifts 
(e.g., Walsh functions ) are applied to each of the input antennas in order to reduce the response of the instrument to 
spurious signals fe.g.. lThompson et al.|[200lh . The phase shift patterns we require in the adding interferometer must obey 
a more stringent orthogonality requirement: rather than merely demanding orthogonality of all of the input phase shifts 
(i.e., demanding that the <fij be orthogonal), we require that the phase shifts associated with all visibilities (i.e., all 3>jfc) 
be orthogonal. 



l-TjTe-i&it-MSt (9) 



3. Method for finding phase shifts 

Let us suppose that the number N of horns is fixed, as is the number P of possible phase shift values. We wish to find 
the shortest sequence of time steps (that is, the minimum M) that satisfies our orthogonality criterion. Alternatively, 
given M, P, we can ask for the maximum number of horns that can be accommodated. 
We will introduce the following shorthand notation for the possible phase factors: 

\p]=e i2wp/P , P = 0,1,2,...P-1. (10) 

For purposes of illustration, we will consider the case P = 4 in this section, so that the four possible phase shift values 
are 

[0123] = (l,i,-l,-»). (11) 

The method we outline generalizes to other values of P. 

Let the number of time steps M be a power of 4: M — 4 M for some positive integer fi. We can define a set of /i 
mutually orthogonal M-dimensional vectors as follows: the vector a.\ is obtained by stepping through the four possible 
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phase values as slowly as possible - that is, it consists of M/4 repetitions of [0], followed by M/4 repetitions of [1], etc. 
Each subsequent vector a.j cycles through the possible phases four times faster until the last one cc M , which consists of 
M/4 repetitions of the sequence [0123]. To be explicit, here is the case /i = 3: 

ax = [0000000000000000111111111111111122222222222222223333333333333333] (12) 

ol 2 = [0000111122223333000011112222333300001111222233330000111122223333] (13) 

a 3 = [0123012301230123012301230123012301230123012301230123012301230123] (14) 

Here is an alternative description of the construction of these vectors: the fcth element of the vector ctj is the jth-most- 
significant digit in the base-4 expression for (k — 1). In the general case with M ~ P M , a.\ steps through the P values as 
slowly as possible and each subsequent ctj cycles P times faster. 
We now define 

■ ■ ■ , h,ji) = ccj? ■ ■ ■ aj'aj 1 (15) 

for integers , . . . ji > ji between and 3. Here multiplication and exponentiation are performed clcmcntwise in each vector. 
Since [p] is shorthand for e lp7r / 2 , multiplication corresponds to addition modulo 4 on the values in square brackets. For 
instance, in the case \i = 2, 

(2,l) = o^ai = [0123012301230123] 2 [0000111122223333] 

= [0202020202020202] [0000111122223333] = [0202131320203131]. (16) 

It is straightforward to check that the vectors (j^, . . . , j\) are all mutually orthogonal. Since there are 4^ distinct 
vectors, they are a maximal set of orthogonal vectors. We can therefore search among this set for the optimal set of N 
phase shift patterns to apply to our input horns. 

As an example, consider the case /i = 2, that is, let the number of time steps be M = 4 2 = 16. We will determine 
the maximum value of N that can be accommodated. We proceed by assigning phase shift sequences to the horns one 
at a time. Without loss of generality, we can assume that the first horn has no phase shift at all (since any phase shift 
sequence can be subtracted from all inputs without altering the solution): 

00 = (0, 0) = OL%a\ = a Q = [0000000000000000]. (17) 
Here <pj refers to the M -dimensional vector (e 1 ^ 1 , e 1 ^ 2 , . . . , e** > -> M ). We can accommodate two more inputs by choosing 

01 = (0, 1) = a° 2 a\ =ai = [0000111122223333], (18) 
4>2 = (1, 0) = a^aj = a 2 = [0123012301230123]. (19) 
For these three input horns, we have six distinct baselines, with masks [equation ((SJ)] 

*oi = (0,3), $02 = (3,0), * 12 = (3,1), (20) 

*io = (0,l>, * 20 = (1,0), * 21 = (1,3). (21) 

These are obtained by subtracting the values in angle brackets for the two horns modulo 4. for instance, $12 = (0, 1) — 
(1, 0) = (—1, 1) = (3, 1). These masks are all distinct, and hence mutually orthogonal, and furthermore are all orthogonal 
to the vector (0,0), which is sensitive to the total power. 

This construction shows that we can accommodate three horns with a sequence of 16 time steps. We next ask whether 
it is possible to accommodate a fourth vector </>3 in such a way that the new masks $03 > &13 > etc. are independent of the 
ones we have already found. A search of the 16 — 7 = 9 candidates reveals an affirmative answer: <p 3 = (3, 3) = 0^20:1 
works. 

The value N — 4 is the maximum that can be achieved for the case of M = 4 2 time steps, as is clear from a counting 
argument: N — 5 horns would require at least M = N(N — 1) + 1 = 21 steps. 

Table [1] shows the maximum number of horns that can be accommodated for various values of M. These were found 
by recursively searching the space of possible phase shifts in the manner described above. The last column shows the 
maximum value that would be possible according to the simple counting argument that the number of time steps must 
exceed the number of baselines. We have repeated this analysis for the case P — 2, where the phase shifters are capable 
of only and 180° shifts, and found very similar results for the relationship between M and N. 

As noted in the introduction, extremely large values of M are impractical. This is one reason that a bolometric 
interferometer with a large number of horns should surely be designed with a high degree of symmetry, so that there are 
many equivalent baselines that can be read out coherently. (The other reason is the signal-to-noise advantage.) 

In summary, this section has presented a procedure for selecting <fij that yields fully orthogonal masks. This means 
that the result of applying the mask for a given baseline will only be sensitive to the signal from the desired baseline, or 
cquivalcntly that the reconstruction of all visibilities is accomplished with minimal noise. 
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M 


iVmax (actual) 


iV max (counting) 


4 


2 


2 


4* = 16 


4 


4 


4 a = 64 


8 


8 


4 4 = 256 


15 


16 


4 b = 1024 


24 


32 


4 b = 4096 


40 


64 



Table 1. The number of horns N that can be accommodated with a given number of time steps M. We assume P = 4 
distinct phase shift values. The second column shows the maximum number that can be accommodated, while the third 
column shows the number found by the simple counting argument N{N — 1) + 1 < M . 



4. Square Array 

We now consider the case where the input horns are arranged in a square array with iV s id. e horns on a side. In this case, 
many different baselines (i.e., pairs of horns) sample the same visibility. We wish to apply identical phase shifts to such 
equivalent baselines, so that a single mask reads out their sum. Naturally, we also require that inequivalent baselines 
have orthogonal masks. This is the case considered in detail by Charlassier et al. (C08). Our method parallels theirs in 
many respects but refines it in some ways. 

Following the notation of C08, we parameterize the position of horns in the array in units of the minimum horn 
separation as a vector dj = (lj, rrij). Here lj, rrij are integers running from to iV S id c ~ 1, labeling the position of the horn 
in along the x and y directions. The index j runs from to N 2 idc — 1 according to j = lj + N S id e mj. We can construct 
a set of phase shifts for all horns that satisfy the desired criteria using the basis described in the previous section with 
[i = 2. We let the phase shift sequence for horn j be 

4>j = Qi>mj). (22) 
Below we have explicitly written out the modulations for each horn in a 6 x 6 array. 

(0,0) (1,0) (2,0) (3,0) (4,0) (5,0) 

(0,1) (1,1) (2,1) (3,1) (4,1) (5,1) 

(0,2) (1,2) (2,2) (3,2) (4,2) (5,2) 

(0,3) (1,3) (2,3) (3,3) (4,3) (5,3) 

(0,4) (1,4) (2,4) (3,4) (4,4) (5,4) 

(0,5) (1,5) (2,5) (3,5) (4,5) (5,5) 



In this case, the mask for the visibility corresponding to any pair of horns is simply (AZ, Am). This means that all 
pairs with the same relative spacing get the same mask. Furthermore, as long as the number of phase shift steps P is 
large enough, all inequivalent visibilities correspond to orthogonal masks as desired. The minimum value of P is set by 
the fact that phases are only defined modulo P. Since Al, Am can range from — (iV 8 ide — 1) to N S id e — 1 we need at least 
P = 2(iV s idc — 1) + 1 distinct phase shifts. (As we will see in the next section, it may be desirable for P to be a multiple 
of 3, in which case we simply round up to the nearest such value.) If P is smaller than this, then distinct visibilities will 
be mapped onto the same phase shift sequence. For the above case, for example, we require P > 11. If we tried a smaller 
value, say P = 10, then the visibility corresponding to horns (0,0) and (5,1), for example, would get the same phase shift 
sequence as (5,0) and (0,1), namely (5, 1) = (—5, 1). 

It is instructive to compare this scheme with the very similar one of C08. In both methods, the phase shift sequence for 
horn (I, m) is expressed in the form Ih + mv for two basis shift patterns h, v. In order to achieve the desired orthogonality 
properties, Charlassier et al. choose h, v to be independent random vectors of phase shifts. The randomness ensures 
approximate orthogonality, up to errors of order M~ x l 2 , where M is the length of the phase shift sequence. In contrast, 
we choose h = (1,0) and v = (0, 1). This results in strict orthogonality, as opposed to approximate orthogonality. 

The number P of distinct phase shift values required is essentially the same in the two methods. As in our method, 
C08 found that P ~ 27V s id e was required in order to produce orthogonal phase shifts using random basis vectors. 

Using either method, the length M of the modulation sequences is greatly reduced compared to the case of inequiv- 
alent baselines. The number of required phase shifts is M = P 2 = (2N S ide — l) 2 when redundant baselines are tagged 
equivalently. If we instead used the methods of the previous section, we would require M > N 2 idc (N 2 idc — 1) + 1 w A r s 4 idc . 
For the 6x6 array denoted above, this is the difference between a 121-step sequence and a 1261-step sequence. Even 
more important is the signal-to-noise benefit of coadding equivalent baselines. 

Although we have described this procedure as applying to a square array, it is in fact more general. It applies whenever 
the horn positions can be expressed as integer multiples of any two basis vectors, even if the two are not orthogonal, or 
in other words, to any parallelogram-shaped array. Furthermore, it can be applied to any subset of such a parallelogram- 
shaped array, since we can simply ignore the parts of the array with no horns in them. In particular, this means that the 
method can be applied to a hexagonal close-packed array of horns, as shown in Fig. [TJ 
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Fig. 1. A hexagonal array of 19 horns can be seen as a subset of a 5 x 5 parallelogram-shaped array. The phase shifting 
scheme described in Sect. [4] with A^idc = 5 can be applied to this array. 



5. Two Polarizations 

Thus far, for simplicity we have been considering only one polarization state of the incoming radiation field. We now 
imagine that two orthogonal linear polarizations (x, y) are measured at each horn. In this case, we can in principle 
recover visibilities for all four Stokes parameters I, Q, U, V. To be specific, if visibilities V xx , V xy , V yx , V yy are measured 
for a particular baseline, then 



Vj = V XX + Vyy, (23) 

V Q = V XX - Vyy, (24) 

V V = V X y + Vy X , (25) 

V v = -i(V xy ~ V yx ). (26) 



However, as C08 have pointed out, it is impossible to recover all visibilities while taking full advantage of the noise 
reduction resulting from coadding redundant baselines. C08 describe two schemes for recovering some of the Stokes 
parameters with full accuracy, one of which (mode 2 of C08) involves measuring the visibilities for Stokes /, U, V but not 
Q. Stokes Q can then be measured by rotating the instrument 45°. In this section, we show how to implement this mode 
of operation using our phase shifting scheme. 

Aside from the phase shifting scheme, there is another reason for adopting an observing scheme in which Stokes Q is 
measured only by rotating the instrument. As eq. (j2"4"|) indicates, the visibility for Stokes Q is obtained by subtracting 
two measured visibilities, each of which contains a contribution proportional to the much larger Stokes parameter /. As 
a result, this visibility is likely to be subject to much larger errors than the other linear polarization (Stokes U). 

As in the previous section, we assume an A s ;d G x N s a c array of horns, but now we introduce an orthomode transducer 
for each horn, doubling the number of signals to be interfered. We can represent each of these 2A 2 idc signals with a triple 
of labels (lj,rrij,nj) where (lj,rrij) label the position of the horn as in the previous section, and rij =0,1 labels the 
polarization state. In the previous section, we identified the horn (lj,rrij) with a phase shift sequence (lj,mj). In the 
present case, we can define phase shift sequences similarly in terms of the triple {lj, rrij, rij), each of which represents a a 
sequence of 3P 2 time steps, where P > 2A S jd — 1 as in the previous section. For one of the two polarization states, we 
the phase shift sequences are simply three repetitions of the sequences we found previously: 

(I, m, 0} = ((/, m), (I, m), (I, m)) . (27) 

For the other polarization state, we apply a slow three-phase modulation to this sequence: the first P 2 steps are unchanged, 
the next P 2 steps are multiplied by e 27 ™/ 3 , and the final block is multiplied by e 47 ™/ 3 : 

(Z,m,l>= ((?,m),e 27rl/3 (Z,m),e 47r4/3 (/,m)) . (28) 

Note that this scheme is most natural to apply when P is a multiple of 3 so that for every phase state p there is another 
whose phase is p + 2ir/3. Otherwise, the set of phase shifts involved in the sequences {l,m, 1) will be larger than that 
involved in (l,m, 0). In implementing this scheme, one would surely round P up to the nearest multiple of 3. 

It is straightforward to check that all equivalent baselines have identical phase shift sequences as desired. All pairs 
that interfere x and y polarization have independent, orthogonal phase shift sequences, allowing optimal reconstruction 
of Stokes U, V visibilities [eqs. (j25"]) . ([26"l) ]. Those that interfere x and x have identical sequences to those that interfere y 
and y. Applying these phase shift masks therefore allows recovery of the sum of these visiblities, which is Vj . 

As in the previous section, this method is similar to that of C08, except that our method imposes strict orthogonality 
on distinct baselines, as opposed to relying on the approximate statistical orthogonality that results from choosing random 
phase shift sequences. 

As an example, consider a square array with N s - l( io = 8. The number of different phase shift values must satisfy 
P > 2A S ido — 1 = 15. The length of the phase shift sequences is M = 3P 2 = 675. The shortest sequence of phase shifts we 
could possibly hope for would have M equal to the number of unknowns we are trying to solve for. In this arrangement, 
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there are 112 incquivalent baselines, each of which has three complex visibilities that are measured, and in addition the 
total power in /, Q, U are measured, resulting in a total of 6 x 112 + 3 = 675 unknowns. Our phase shift sequence is 
therefore as short as possible. For comparison, according to Fig. 4 of C08, the optimal noise levels in the C08 scheme are 
obtained only when the phase shift sequence is ~ 3 times the minimum length. In these comparisons we only consider 
the length of mode 2 in C08. It should be noted that when 27V s id c — 1 is divisible by 3 that we recover visibilities with 
maximum efficiency. When this is not the case the ratio of our length to the minimum approaches 1 for large iV s ide- For 
arrays of reasonable size (8 x 8 or larger) the maximum ratio is 1.22 and occurs for a 10 x 10 array. 

6. Conclusions 

Optimal recovery of visibilities in a bolometric interferometer depends on the proper choice of phase shift sequences. We 
have laid out a method for finding such sequences that lead to fully orthogonal masks for all visibilities and introduced 
a compact notation for describing such phase shift sequences. 

In the case of an array with a regular lattice structure, equivalent baselines can be read out simultaneously, reducing 
the length of the required phase shift sequence and improving the signal-to-noisc ratio. This method refines that of C08. 
The method applies to arrays that are based on replication of any parallelogram-shaped fundamental cell, including for 
example hexagonal arrays. 

For the case of rectangular arrays, the method described herein is very similar to that of C08, although our method 
imposes strict orthogonality on the masks for distinct baselines, rather than relying on approximate orthogonality resulting 
from random sequences. As a result, our method leads to optimal recovery of visibilities for Stokes' /, U, V, with shorter 
time sequences than that of C08. 

The ability to shorten the sequence of phase shifts is likely to be important in instrument design, because it reduces 
the degree to which 1// noise must be controlled. For example, suppose that we can shift phase states at a rate of 
one state per 10 ms (either because of the design of the phase shifters or the bolometer time constants). As we saw in 
the previous section, an 8 x 8 array requires ~ 1000 phase shifts, which would take 10 seconds. We therefore require 
the 1// noise knee to be below ~ 0.1 Hz. An alternative scheme involving a longer phase shift sequence would require 
correspondingly tighter control of the 1// knee. 
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Appendix A: Proof of minimum-noise condition 

In this section we provide a proof of the assertion that orthogonal phase shift patterns minimize the noise in the recovered 
visibilities. 

Assume that the visibilities are arranged in an A^s-dimensional vector V, and the observed signals are arranged in 
an M-dimensional vector S. The two are related by an M x iV v i s matrix A: 



All entries of A are complex numbers with absolute value 1. We assume that M > N v i s and that the matrix A has 
maximal rank, so that it is possible to solve for the unknown visibilities. 

Assuming that the signals are contaminated with white noise with variance er 2 , the optimal reconstruction of the 
visibilities is the least-squares vector 



S = AV 



(A.l) 



V = (A^A^A^S. 



(A.2) 
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The noise covariance matrix for V is 

Af = a 2 (A^A)~ 1 . (A.3) 

The noise in the jth recovered visibility has variance Afjj . We wish to show that this noise is minimized when the matrix 
A has orthogonal columns. 

The diagonal elements of the inverse noise matrix arc 

M 

(Af-^jj = a-'(AU) i i = a- 2 ]T A* m3 A mj = . (A.4) 



m— 1 



We can therefore write 



7V- 1 = ^(J + ©), (A.5) 

where I is the identity matrix and V is a hermitian matrix with zeroes along the diagonal. 
The noise covariance matrix is 

2 2 2 

N=— (1-V + V 2 -V z + ...) = — \1-V + V(l-V + V 2 - .. .)£>] = —(1 — 7) + VMV). (A.6) 
Since V has no diagonal elements, an arbitrary diagonal element of the noise covariance matrix is 

*ii = ^[1 + WV)^] = ^[1 + vWv], (A.7) 
where Vk = Dkj- Since M is a positive definite matrix, we conclude that 

MjJ > ^- (A.8) 

That is, the minimum noise variance achievable on any one visibility is a 2 /M . This value is achieved when the matrix A 
is column orthogonal, since in this case A^A — (M/a 2 )l and J\f = (a 2 /M)I. 



